Dynamics of quantum-classical hybrid system: effect of matter-wave pressure 
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Radiation pressure affects the kinetics of a system exposed to the radiation and it constitutes 
the basis of laser cooling. In this paper, we study matter-wave pressure through examining the 
dynamics of a quantum-classical hybrid system. The quantum and classical subsystem have no 
explicit coupling to each other, but affect mutually via a changing boundary condition. Two systems, 
i.e., an atom and a Bose-Einstein condensate(BEC), are considered as the quantum subsystems, 
while an oscillating wall is taken as the classical subsystem. We show that the classical subsystem 
would experience a force proportional to Q -3 from the quantum atom, whereas it acquires an 
additional force proportional to Q~ 2 from the BEC due to the atom-atom interaction in the BEC. 
These forces can be understood as the matter-wave pressure. 

PACS numbers: 73.40.Gk, 03.65.Ud, 42.50.Pq 



I. INTRODUCTION 

It is well known that electromagnetic radiation exerts 
a pressure upon any surface exposed to it, the pressure 
was deduced theoretically by James Clerk Maxwell and 
Adolfo Bartoli more than a century ago, and proven ex- 
perimentally by Lebedev, Ernest Fox Nichols and Gor- 
don Ferric Hullin in the last century. Recently, this re- 
search field becomes activeQ, H, i, i, H, B 0, B again 
due to the breakthrough in nanofabrication and in ul- 
tracold science, this together with the coupling of co- 
herent optical system to micromechanical devices, has 
opened up the exciting new field of research, cavity 
optomechanics[I3, [U El El El • 

Quantum mechanics tells us that all matter exhibits 
both wave-like and particle-like properties, called wave- 
particle duality, the wavelength is inversely proportional 
to the momentum of a particle and the frequency is di- 
rectly proportional to the particle energy. These facts 
give rise to a question: Can the matter wave exert a 
pressure (matter-wave pressure) upon a surface exposed 
to it, like the electromagnetic radiation does? To answer 
this question, it is good to study a coupled quantum- 
classical hybrid system, and examine how the quantum 
subsystem affects the kinetics of the classical subsystem. 

On the other hand, classical-quantum hybrid system 
on its own is an interesting topic to study [IH Il6| . By 
classical-quantum hybrid system we mean a composite 
system consisting of a quantum and a classical subsys- 
tem. For a closed hybrid system, the quantum subsystem 
may be treated classically for specific issues addressed 
[l7l[T§|. However, this approach in general is inadequate 
owing to the loss of quantum features. In Ref.[l5[, the 
authors present a general framework for exact treatment 
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of such a hybrid system. When the quantum subsys- 
tem is dynamically fast and the classical subsystem is 
slow, a vector potential arises. This vector potential, on 
one hand, gives rise to the familiar Berry phase in the 
fast quantum dynamics, on the other hand, it yields a 
Lorentz-like force in the slow classical dynamics. In the 
formalism, the Hamiltonian of the fast quantum subsys- 
tem depends explicitly on the freedom of the classical 
subsystem. In contrast, here we study the dynamics of 
quantum-classical hybrid system without any explicit in- 
teraction, instead an boundary condition will be consid- 
ered that connects the quantum and the classical subsys- 
tem. An atom and a BEC are taken as the quantum sub- 
system, both of them move in an one-dimensional square 
well with an oscillating wall as their boundary. The dif- 
ference between these two quantum subsystems is that 
the BEC has atom-atom interaction, whereas the atom 
would not have. This setting is exactly a model we need 
to discuss the problem of matter-wave pressure. 

The paper is organized as follows. After a brief intro- 
duction to the formalism for classical-quantum hybrid 
system in Sec. II, we study the kinetics of the classical 
subsystem (i.e., the oscillating wall) under the effect of a 
quantum atom in Sec. III. Treating the atom as a classi- 
cal particle, the kinetics of the wall is also examined in 
this section. In Sec. IV, we investigate the kinetics of the 
wall exposed to a Bose-Einstein condensate. The atom- 
atom coupling in the BEC will make the force acting on 
the wall different from the case in Sec. III. Finally we 
conclude our results in Sec.V. 



II. GENERAL FORMALISM 

In this section, we present a brief introduction to 
the formalism for quantum-classical hybrid system. For 
more detail, we refer the reader to Ref . [l5j . Consider a 
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quantum-classical hybrid system described by 

if =<*l#i(Q)l*}+ff 2 (P,Q), (l) 

where Hi is the Hamiltonian for the quantum subsystem 
and I*) = (ipi, ip 2 , ^Pn) T denotes its quantum state. 
H 2 (P,Q) represents the classical subsystem and P, Q 
stand for its coordinates and momenta, respectively. It 
has been shown [i~9H20j that a quantum system possesses 
mathematically a canonical classical Hamiltonian struc- 
ture. This can be understood by defining Xi — \fift}$>i, 
iji = \fiftM>*. With this definition, one can rewrite the 
dynamics for the quantum system as 



dxi 
lit 

dyi 
dt 



dHj_ 

_dH 1 
dxi ' 



(2) 



where H[ = H^x^Q) = #i(*,**,Q) = 
(\I , |ffi(Q)| 1 4'). Thus the quantum system can be refor- 
mulated in the language of classical theory Eq. ((2|) . 

The quantum state can be expanded in terms of 
instantaneous eigenstates of the Hamiltonian Hi, 



N 

E< 

n=l 



i(Q)>, 



(3) 



where H 1 (Q}J^„(Q)) = E n (Q)\<j> n (Q)) . By the formal- 
ism in Ref. [15( , we can obtain a vector potential that the 
classical subsystem feels, 



A' 



8 



A = Ji n A ni A n = i(<^ n | — |^ n ) 



(4) 



where Ii n — h\c n \ 2 . This leads to a dynamical equation 
for the classical subsystem, 



MQ 



dJf°i 



dQ aq 



dv 2 . 

T^T + Q x 



(5) 



with tfPx = (Ii, Q) = J2 n E n (Q)I ln /h, V 2 is a poten- 
tial and 33 = V x A = J2 n hn^ x A„ is a gauge field 
like magnetic field. 



III. ATOM AS THE QUANTUM SUBSYSTEM 

Consider a bipartite hybrid system which consists of 
a single atom in an infinitely deep well and a moving 
wall, the moving wall acts as a boundary for the atom 
(see FiglT|), and the whole system is restricted to move 
in one-dimension. The Hamiltonian of such a system can 
be described by 



H = Hi + H 2 



(6) 



M , Q CO 




FIG. 1: Schematic diagram of the quantum-classical hybrid 
system. 



where Hi denotes the Hamiltonian of the atom trapped 
in the infinitely deep well, H 2 stands for the Hamilto- 
nian of the moving wall which is considered as a classi- 
cal harmonic oscillator. Here Hi = p 2 /2m, and H 2 = 
■J^j + \Mtu 2 (Q — Qo) 2 - P an d P denote the momen- 
tum along the x axis for the quantum atom and classical 
wall respectively; Note p is an operator while P is a c- 
number here; m and M stand for the mass of the atom 
and the oscillator respectively; V 2 = \Mlo 2 {G> — Qo) 2 
is the potential energy of the harmonic oscillator; Q is 
the coordinate of the moving wall and Qo is its equilib- 
rium position; ui denotes the vibration frequency of the 
classical subsystem. In this hybrid system, we consider 
a situation that there is no explicit interaction between 
the quantum atom and the moving wall. The moving wall 
affects the dynamics of the quantum atom by changing 
its boundary condition. Interesting features arise in this 
case as we will show below. 

In our model the classical subsystem only changes the 
boundary condition of the quantum subsystem, which 
would reflect in the instantaneous eigenstate and the cor- 
responding eigenvalue of the quantum subsystem given 
below 



\<Pn{Q)) = / (j>n(q,Q)\q)dq, 

Jo 

a, l n\ I** ■ n7Tq 

<Pn{q,Q) = y Q sm ~Q~ > 



E n (Q) 



h 2 n 2 n 2 
2mQ 2 ' 



(7) 



where \q) denotes the eigenstate of coordinate q. Then 
the Hamiltonian of the quantum subsystem can be 
rewritten as 



q) = E n {Q)hn/n = J2 



n 2 ir 2 h 2 
2mQ 2 



(8) 



With these results, the vector potential A is then 
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FIG. 2: The coordinate of the classical subsystem Q (in units 
of nm) as a function of time (in units of 10 _7 s) with B = 0.01, 
lu = 1 (in units of 2nx 10 7 Hz), Qo = 1, the initial condition is 
Q(0) = 1.1 and Q(0) = 0. In the red solid curve, we consider 
the effect of the quantum subsystem, while in the blue dash 
curve we plot the free evolution of the oscillator. 

■A f Q d 

= 2_.iH\c n \ 2 / <t>* n { ( l,Q)-^ ( t>n{q-,Q)dq 

= 0, (9) 

leading to an observation that only scalar potential (from 
the quantum subsystem) affects the dynamics of the clas- 
sical subsystem. It is worth noting that this conclusion 
depends on the model, i.e., the magnetic-like gauge fields 
are not zero in general. 

Substituting these equations together with the poten- 
tial energy V% into Eq.©, we get the dynamical equa- 
tion for the classical subsystem of the linear system (one- 
dimensional) 

Q = -^-u 2 (Q-Q ), (10) 

2 2^2 

where B = J2 n n mM • This equation includes the effect 
of the quantum subsystem and describes the kinetics of 
the classical subsystem. The first term in Eq. ([T0|) rep- 
resents a force acting on the wall from the atom. We 
explain this force as a consequence of matter-wave pres- 
sure for the following reasons. (1) Suppose that the pop- 
ulation of the atom in level n remains unchanged, the 
total energy of the atom is E\ — £ n "2^ l c ri| 2 , where 
\c n \ is the probability of the atom in level n. Consider 
a small vibration AQ of the wall, the work done by the 
atom is F ■ AQ, and we have F = — ^q, this exactly 
leads to the first term in Eq. (fT0f . (2) The wave- function 
of the atom is |$) = £„ c„|<^„(Q)), with the energy 
Ei and momentum p = J2 n \c n \ 2 Pn = J2n M 2 ^- the 
probability current is J = £ n J ™ = J2 n \ c ™\ 2 ^q- Tnis 
means that the atom exerts a pressure force given by 
F = En 2p„J„ = £„ \cn\ 2 ^^ on the wall. 

In FigfJl we plot the coordinate of the classical sub- 
system Q as a function of time. It is clear that Q is a 
oscillating function of time. The blue-dash line denotes 
a free harmonic oscillator and the red-solid line is for 
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FIG. 3: (Color online) The time evolution for the coordinate 
of the classical subsystem Q under different initial amplitudes 
(Q(0) — Qo) f° r the parameters of Q(0) = 1.1 (red solid curve) 
and Q(0) — 1.8 (blue dash curve). The other parameters 
are chosen as B — 0.1, u = 1, Q(0) = and Qo — 1. The 
frequency, time and length were rescaled as the same as Fig[5] 

the case with the effect of the quantum subsystem. Here 
the parameters are chosen as B = 0.01, u> = 1, Qo = 1 
Q(0) = 1.1 and Q(0) = 0, where Q(0) and Q(0) denote 
the initial maximal displacement from the equilibrium 
and the initial velocity of the moving wall, respectively. 
Comparing with the free harmonic oscillator (the blue 
dash line in Figj2j), we find that when we consider the ef- 
fect of quantum subsystem the equilibrium point moves 
to the right side. This is a consequence of the matter- 
wave pressure. In Fig[3l we present the time evolution 
of the coordinate of the classical subsystem Q under dif- 
ferent initial amplitudes (Q(0) — Qo)- For the red solid 
curve, we choose Q{0) = 1.1, while for the blue dash 
curve, Q(0) — 1.8. The other parameters are B = 0.1, 
uj = 1, and Qo = 1- We find that the effect of the 
quantum subsystem increases with Q(0) — Qo- The curve 
becomes sharper when the coordinate moves toward the 
origin. This feature becomes evident for the case of large 
initial amplitude. Next we show the coordinate Q as a 
function of time with different initial equilibrium position 
Qo in FigHJ We find that for the same initial amplitude, 
the larger initial equilibrium position Qo is, the smaller 
the effect of quantum subsystem on the classical subsys- 
tem. 

These observations can be understood by examin- 
ing the dynamical equation Eq. (|10p . In this equation 
MQ denotes the resultant force of the moving wall and 
Muj 2 (Q — Qo) is the elastic spring force which keeps the 
moving wall in harmonic oscillating. Quantum subsys- 
tem brings in a rightward force (^§-, with Bo = BM) 
which is inversely proportional to the cube of the coor- 
dinate Q of the moving wall. This is the reason why the 
equilibrium position moves more to the right side than 
the free harmonic oscillator. The quantum subsystem in- 
duced force results from the matter-wave pressure, hence 
by examining the kinetics of the classical subsystem, we 
may recognize some features of the matter- wave pressure. 
It is easy to see that this force decreases rapidly with the 
increasing of Q. This is the reason why the effect of the 
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FIG. 4: (Color online) The time evolution of the coordinate 
Q for different initial equilibrium positions Qo but the same 
initial amplitude (Q(0) — Qo = 0.1). Here the parameters are 
chosen as B — u = 0.1, Q(0) = 0. The red solid curve is 
plotted for Qo = 3, while for the blue dash curve, Qo = 1. 
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FIG. 5: (Color online) In contrast to Figj4j in which that atom 
was treated as a quantum system, in this figure we show the 
coordinate Q (thick green) and q (thin blue) as a function of 
time. Both the atom and the wall are considered as classical 
objects. This plot is for different initial equilibrium positions 
Qo, but the same initial amplitude (Q(0) — Qo = 0.1), (a) 
Qo = 3, (b) Qo = 1. The other parameters chosen are Mu> 2 = 
60, M/m — 1000, the initial velocity of the atom and the wall 
are 25 and 0, respectively. The atom was assumed initially 
at the origin. In this plot, time is in units of 10 _7 s, length in 
units of nm. 



quantum subsystem become smaller when the width of 
the well gets larger. In other words, this force changes 
in a wide range when the amplitude of the oscillation (of 
the wall) is large. 

It is interesting to compare these observations with the 
results given by classical mechanics. To this end, we sim- 
ulate numerically the kinetics of the hybrid system, treat- 
ing both the atom and the wall classically. The Hamilto- 
nian that governs the kinetics of such a classical system 
is the same as Eq.(|6]), but p in Hi is not an operator 
now. The atom is bouncing back and forth and colliding 
with the wall. We assume that the collision is elastic such 
that both energy and momentum are conserved. The nu- 
merical results are presented in FigfSl where we plot the 
coordinates of the atom (q) and wall (Q) as a function of 
time. Similar to the quantum case, the smaller the well is 
(or, the smaller the Qo), the lager the effect of the atom 



on the wall. The difference is that a quantum atom exert 
a force proportional to 1/Q 3 on the wall, while a classi- 
cal atom provides a force on average proportional to 1/Q 
according to the law of conservation of momentum. As a 
consequence, the kinetics of the wall behaves differently 
for Q -> 0. 



IV. BOSE- EINSTEIN CONDENSATE AS THE 
QUANTUM SUBSYSTEM 

In this section we will take a BEC as the quantum sub- 
system to study the kinetics of the classical subsystem. 
The difference between the atom and the BEC that will 
manifest in our study is the atom-atom interaction in the 
BEC. We will show in the following that this atom-atom 
interaction results in an additional force proportional to 
Q~ 2 to the classical subsystem. For simplicity, we con- 
sider a BEC in one-dimensional square well with an mov- 
ing wall as its boundary. The setting is the same as in 
Figffl but the atom is replaced with a BEC. The Hamil- 
tonian which describes such a system can be written as, 



H 



ij*(x)H 1 (x, Q)^(x)d 3 x + H 2 {P, Q), (11) 



where Hi denotes the Hamiltonian for the BEC, and 
ip{x) is the wave function of the BEC. H2{P,Q) is the 
Hamiltonian of the classical moving wall and P, Q are 
its coordinate and momentum, respectively. The effective 
Hamiltonian of a BEC in a potential V(x) takes, 



ti 2 d 2 



g\ip(x,t)\ 2 , 

where the potential V(x) in our model is, 
V(x) = 



0<x<Q 
f-oo x < or x > Q 



(12) 



(13) 



Here x stands for the coordinate of the BEC and the 
atom-atom coupling constant in BEC is denoted by g. 
The stationary CP equation can be written as 

ft 2 d 2 \ 

-^-7Ti + + 9\i>( x > t)\ 2 V'Oz, t) = wl>(x, t). 
2m ox z ) 

(14) 

Here u denotes the chemical potential. In terms of Jacobi 
elliptic functions, the eigenfunction for repulsive interac- 
tiono > and attractive interaction g < can be written 

as nin, 



<p+A x ) 

<p-,j(x) 



6 +J sn(a +J • x + S+j, k +J ), 
&_jcn(a_j • x + S-j, fc_j), 



with 



for g > 0, 
for g < 0, 
(15) 



(16) 
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where sn, en are the Jacobi elliptic functions and 
k±,j C? = 1,2,3,...) are the modular number of the Ja- 
cobi elliptic function, S±,j (j — 1,2,3, ...) are constants 
which will be given below. We use + and — to denote 
the case of g > and g < respectively, and j labels the 
eigenf unctions. 

Taking the boundary conditions p±.j(0) = (f± t j(Q) = 

and the normalization condition J Q \ip±j(x)\ 2 dx = 1 
into account, we obtain, 



2jK(k 



l +,3 



3 = 1,2,3. 



b +.j = 



k +tj K(k +J ) 



[K(k +tj )-E(k +tj )}Q> 

6- d = -K(k- d ), 

2jK(k- ,•) 
a-j= J K Q ' 3) , j = 1,2,3... 



[E[k- d ) + (fc_j - l)Jf(fc_j)]Q 



. (17) 



where K(k± t j) and E(k±j) are the first and the second 
elliptic integrals, respectively. 

We now focus on the kinetic equation of the classical 
subsystem. In our model, it is easy to show that SS — 
V x A = 0, and the classical Hamiltonian for the BEC 
takes 



M 1 = 



Q ,*/ & 9 2 

2m ox 



V(x) 



1 



+t;9V* {%, t)v(x, t)]y(x, t)dx. 



(18) 



Specifically, for g > 
and for g < 

h 2 



2jk +J K 2 {k 



2jK(k^ 



2m Q 2 [K{k + ^) - E(k+j)] J 
gk 2 +J K(k + . J ) 



4j[K(k +J ) - E(k+j)]*Q Jo 



l)sn 2 z — 2kj r jsa A z\dz 



sn zdz, 



2jk-jK 2 (k = j) 



2m Q 2 [{k-^ - l)K(k- tj ) + E(k- tj )] 
+ ij[(k-j-l)K(k-j)+E{k-^Q 



(2j-l)K{k-j) 

[1 - (2k— j + l)sn 2 y + 2k_ d sn 4 y]dy 

-K{k- :j ) 

(2j-l)K(fe_, J ) 

cn A ydy. (g < 0) 



(19) 



(20) 



Here z = 2 ^ K ^+^ X and y — K(k__j)( 2 0- — 1). Substituting these equations together with potential V — jMuj 2 (Q — 
Qo) 2 into Eq.([5]), we obtain a kinetical equation for the classical subsystem. For repulsive interaction, i.e., g > 0, it is 



(21) 



with 



and 



Ci 



2jh 2 k +d K{k +tj f 
Mm[K(k+,j) - E(k + ,j)} Jo 



2jJf(*+,j) 



[(fc+j + l)sn 2 z — 2fc +)J sn 4 z](iz > 0, 



gk 2 K(k. 



+ ,3> 



2jM[K{k, 



E(k 



sn zdz > 0. 



For 3 < 0, 



( 3 = § + ^|-- 2 (Q-^), 



(22) 



where 



C 2 = 



2jh 2 k_ J K(k- J ) 2 



(2j-l)Jf(fc_, j ) 



[1 - (2fc_j + l)sn 2 y + 2k-jsn A y]dy > 0, 
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and 



gk 2 _ tj K(k- d ) / .(2j-i)K(fe_,,) ^ 
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FIG. 6: (Color online) The time evolution of the coordinate 
Q (in units of nm) for different atom-atom interactions in 
the BEC. The parameters and initial conditions are chosen 
as Ci = C 2 = 0.01, w = 1 (in units of 2vr x 10 7 Hz), the 
equilibrium position is Qo = 3, the initial position is Q(0) = 
3.1, and Q(0) = 0. 

It is easy to show that if there is no interaction between 
the atoms in BEC (namely, g = 0), C; and Di (i = 1,2) 

reduce to C x = C 2 = 3 ^m&, and D x = D 2 = 0. This is 
exactly the result given in Sec. Ill, where the atom was 
taken as the quantum subsystem. Observing Eqs. (|21|) 
and (f2"2")l . we find that the classical wall experiences a 
force proportional to Q -3 , which is the same as that 
we discussed in Sec. III. In addition to this force, a force 
(~ Q 2 ) inversely proportional to the square of the co- 
ordinate of the classical subsystem appears. This force is 
due to the interaction between atoms in the BEC, which 
is different from the result we discussed in Sec. III. 

In Fig [5] we plot the coordinate of the classical subsys- 
tem Q as a function of time for different D\ and D 2 . Here 
we choose the parameters as C\ = C 2 = 0.01, u> = 1, 
Qo = 3, Q(0) = 3.1 and Q(0) — 0. In contrast, we 
plot the situation without atom-atom interaction as the 
green solid line. From the figure, we find that the equi- 
librium point of the classical subsystem moves right for 
the case of D\ = 2 and D\ = 0.5, whereas for the case of 
D 2 = —2 and D 2 — —0.5, the equilibrium position moves 
left. These can be understood by analyzing the dynam- 
ical equation Eqs.([2"Tf and ([22]). In Eq. ([2"Tj) MQ is the 
resultant force of the moving wall and Mlo 2 (Q — Qo) de- 
notes the elastic spring force which keeps the moving wall 
in harmonic oscillation. Quantum subsystem brings in a 



rightward force (^r) similar to the case in Sec. III. An- 
other term that is inversely proportional to the square of 
the coordinate(|p) comes from the repulsive interaction 
between atoms in the BEC. It results in a repulsive force 
for the moving wall, too. There together can explain why 
the equilibrium position of the wall moves right with re- 
spect to the case without atom-atom coupling. For simi- 
lar reasons the equilibrium position of the wall moves left 
for attractive atom-atom interaction. 



V. CONCLUSION AND DISCUSSIONS 

The dynamics of quantum-classical hybrid system has 
been studied in this paper. Two quantum subsystems are 
taken to discuss the kinetics of the classical subsystem. 
When the quantum subsystem is an atom, the classical 
subsystem experiences a force proportional to its 

distance Q to the fixed wall, meanwhile the energy of the 
atom has been changed because the classical subsystem 
provides a moving boundary, even if the atom remains 
in the same level in time evolution. When the quantum 
subsystem is a BEC, the BEC would exert an additional 
force proportional to Q~ 2 to the moving wall. This 
force comes from the atom-atom interaction in the BEC, 
hence it can act as a witness of the nonlinearity in BEC. 
With current technology, a SiN membrane of effective 
mass M = 4 x 10~ 13 Kg is possible in laboratory, a 
vibration frequency 27r x 1.3 MHz sets the time scale of 
the dynamics to - 10 _7 s. m - 10 _27 Kg, M ~ 10~ 13 Kg, 
and the eigenstate index n ~ 500 may lead to B ~ 0.1(in 
units of (Hz) 2 m 4 ) in the first example, with which the 
kinetics of the classical system has been changed sharply. 
This estimation is conservative. In fact, the frequency 
of the membrane can be ~ GHz, leading to a time scale 
~ 10 _9 s. To our best knowledge, this is the first time 
to show theoretically the effect of matter-wave pressure, 
though the experimental observation is a challenge task. 
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